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Abstract 


IVo mclhods for tliC computation of qooid undulations usin^ potential 
cucfflcicntp and Tx 1* terrestrial anomaly c!a*a are examined. U was found 
that both methods nive the same final result but that the method suRKCSted by 
Molodenskli allo\\s a more simplified error analysis *han the mefhod used by 
Vincent and Marsh. 

Specific equations were considered for the effect of the mass of the 
atmosphere and a cap dependent zero-order undulation term was derived. 
A!thou<;h a correction to a gravity anomaly for the effect of the atmosphere 
is only about -0. 'T nigal, this correction causes a fairly large undulation 
correction ic.g. 2. 3m with a cap size of 20®) that lias not previously been 
considered. 

The accairacy of a geoid undulation computed by these techniques was 
estimated considering anomaly data errors, potential coefficient errors, and 
truncation (only a finite set of potential coefficients being used ) errors. It 
was found that an optimum cap size of 20® should be used. 

Thr ;coid and its accuracy were computed in the Geos -3 calibration area 
using the GEMG potential coefficients and l®x 1® terrestrial anomaly data. 

The accuracy of t’ne computed geoid is on the order of 2:2 m with respect to 
an unlcnown set of best earth jjarameter constants. This geoid was com- 
pared to that computed by Vincent and Marsh where we found a systematic 
difference of 3.Um, at undulation difference variance of (2.Gni)^, and a 
maximum difference of 12 meters. 
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Preface* 


Under NASA contract NASn-2l!s4, we are to invcstj^^nte the recovery 
of mean {gravity anomalies from altimeter data that is obtained from Ceos -3. 
In developinj; the methods for such recovery it became clear that it could l^e 
helpful to have an external check on the altimeter determined geoid. This 
check could be obtained by computing geoid undulations from a combination of 
satellite and terrestrial gravity material. This latter set of geoid undulations 
could be used to remove systemaMc bias that might occur in the altimeter data 
due to orbit deter mintit ion inaccuracies and errors in the altimeter itself. 

Although work had been done in the computation of detailed geolds in 
the Geos - 3 calibration area tas well as other areas), no comprehensive ana- 
lysis of the complete theoretical and numerical procedures has been carried 
out. Thus, as one of the first steps in our gravity anomaly recovery work we 
have prejiared this report which attempts to define computational procedures 
for a dctaiUxl undulation computation and its accuracy in a precise way. 
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Introduction 


The purpose of this report is to consider two prcxcdures for the computation 
of detailed geoids using potential coefficients and Txl® terrestrial gravity mater- 
ial. In doing this we are interested in the results obtained from the two methods, 
as well as in the error analysis associated with each. The specific computation 
area is the Geos - H calibration area since for this area an accurate compulation 
of the geoid undulation is needed. For our pur|X)ses the calibration area was con- 
sidered to be Ix'tween 10® and 20® north latitude and from 277® to 297® east longi- 
tude. Although the resulting geoid ma n»»t ije the most accurate available geoid, 
in terms of data used for its computation, we intend tlie analysis given to demon- 
strate an accurate proeedi re for the computation of such a geoid. 

Instead of giving extensive references to past work in this area we shall re- 
strict ourselves to references directly applicable to the current discussion. 

Method A for computing detailed geoid undulation has liecn used extensively 
by Vincent and Marsh (197-11 and Marsh and Vincent (1973). A (iclailed accuracy 
analysis of .Method A w.-’s given by Happ (1973). Method H is described in Molo- 
denskii, et.als. (19G2, page MG) and Heiskanen and Moritz (19G7, page 259) and 
has been used by Groten and Hummel (1971). Details of both methods will be gi- 
ven in Subsequent sections. 


2. Details of Method A and D 


2.1 Method A 


In Method A, the geoid is considered to be composed of three com- 
ponents N\, Nj and N, such that this sum yields the undulation: 


( 1 ) 


N • Ni t Nj + 


Six*cifically we have being the undulation imjilied by a given set of jiotcntlal 
coefficients. The comjxitation of sueh undulations has lioen discussed by Happ 
(1971). There are two methods of interest here. The first, Method One, is liased 
on the solulii'ii ol the geodetic boundary value problem. The result is: 
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( 2 ) 


N, 


CM 



) (C^^cosmX • sInmA) (sino). 


» iifl 


whore; 

CM is th'! geocentric gravitational constant; 

r is ihc gt'occnfric liistaiicc t(» the |M>int ot which the undulation is 

being comiJUlod; 

a is the c(|uatoriuI radius; 

y is the normal gravity at the eomputation point; 


arc* the dilferences l)el\vei-n the :ielu;il potential eoelileients and th«)se 
implied by the adopted reference ellipsoid. 


In 1‘rJiellee we have; 


(2) 

(-'io ■ 


(4) 

^' 4,0 ^\,o(3B5) - 

t’4,0(Wf) 


w’ith all other C and S values equal to their observed values. In our eompufalions 
the reference coefficients were computed from, 


(5) - -Jj / /s” 

•(G) C4,o(«tr) = “J4 Z/iT 


where the J coefficients are (Cook, l‘.)5‘J): 

J, - 2/3 (f(l- i/2) - m/2( 1- 2f/7 t lli^/lU)) 

(7) 

J 4 -4/35 f(l - f/2) ^7f (1 - f/2) - 5m (1 - 2 f/7)) 

with: 

m ••= i/a^(l - f)/ (IM 

and: 

u; is the angular velocity of the earth talvcn as 7. 292A1 514C7 10”'' rad/sec. 
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The solution for r is found by iteration in (8) after the Kcuid potential and 
other quantities in (8) are Riven. The undulation, Nj , is found by dlffcreneinR 
r with the eorresixjndinR r of a apevified reference ellipsoid. Specifically wc 
have: 



• The procedure using (8) and (9) is essentially that used by Vlncen* and 
Marsh. The results from this procedure will bo the same as the results olj- 
tained from (2) provided that consistent constants of G.M, a, to, f, and W'q 
are used. In both cases a zero-order undulation of the geold is taken to bo 
zero. (See section 9 where the removal of this restriction is discussed). 

The No component of (1) is cominitcd in this method as follows: 


(10) N, - JJ (V- Ag.)S(ii>)da 

Oc 

where: 

H is a mean earth radius km); 

G is a mean value of gravity (979. 8 gals); 


-. 1 - 


Aff. Is close to the mran free air anomaly and will lx; discussed 

in a subsequent section; 

A;;, is the mean anomaly implied by <he |K)teiitial coefficients used 

in computing as in tl’); 

Oe is a limited ca|> al>out the com|>utatiun point; 

S(^). is the Stokes' fimction. 


In practice the intcjiration in (10) is replaced by a numerical IntepraMon. 
The value of An, is: 



wh'^TC: 

A is the area in which the mean anomaly AKi is Ijcins determined. 

Specifically in the computation to be ;;iven here A will corresjiond to a l“x 1® 
mean anomaly. The value of Ag, is computed from (Itapp, 1‘JG7): 

(12) Ag, ^ y (f-l) f (r*cosniA 

r — j \r/ L* it 

i-2 


+ S. sinmX) P, (sino). 


The actual evaluation of (11) with (12) was carried out by the analytic inte- 
gration of the cosmX, sin mX terms and the luimencal integration of the 
term. 

The N-, term of (1) is formally given by: 

(13) N, = — 

0~0c 

where: 

O-Oc represents the remaining glolial t ;m> not Int liKk 1 in a,. The cap o,. 

is chosen in such a way tliat N, c.'in lie negle* f» il. A l ap st/.e discussion will lie 
given in section seven. 
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2.2 Mcfh(xl n 

This method has hcen usetl hy (Irotcn and Hummel (1974), for detailed geoid 
eom|>u(ntioD.s. Here \vc write: 


N » nJ ♦ Na + Nt 


where the primes have been used to dlstinglush these values from , N,, and 
N., given In (1). We have: 

(15) N; = ^ Q£(tt.)Ag; 


where: 

is the Volodenskli Iruneution function (lleiBkanen and Marit/., 1957, 

p.2C0) and is given a»: 


Q^(Vo) = 


S(cos^) P (cos V'') sin bdC 


where: 

Oo Ms the cap size of gravity data to be included in Nj and P; are the I.e- 

gendre polynomials. In the computation to l>e given hero wo have use a program 
supplied by M. K. Paul based on his accurate and fast algorillim (Paul, 1973). 

The value of is the £’th degree comjwnent of the gravity anomaly implied 
by a set of jxjtcntial coefficients and is given by: 




4 s siniiiX)!’ (sinC). 
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W'e Ivive for li, : 


’“rrTr 


where the Oe cap has a radius C’o- anomaly Ajf' is the same as in (10). 
The n!^ term can he written as; 


N.' ^ 


This term represents the information alxjve dejfree £,,, thTit has not been 
included in Nj. 


.*).• The Effect of the .Ma.ss of the Atmosphere 


hi dc.eloj)inn the Stokes' equation for comjiutinn the dlslurbin}; potential or 
geoid und'ijati«'n it was assumed that there were no masses extcrmil to the geoid. 
In fact, there are f«)|K)graphic masses and atmospheric masses.. A consideration 
of the topographic masses can be made by alternate solution of the boundary va- 
lue problem. It turns out that the corrections to the results of Stokes' equation 
are small e.xcept in mountainous areas. Speeifieally, the corrections are zero 
in ocean area.s, sueh as bein.; considered here. In any event, procedures for 
handling the topography in this type of computation are discussed In Mortlz (l'J71, 
1‘J75) and are considered negligible U>r this pa|M.*r. 

In considcrbig the atmosphere, we first define the |H)tential of the geoid to be 
the sum of the gravity jKitcntial of the solid earth and oceans plus the potential of 
the atmosj)l)erc. (’onscquently, the computati<»n of tlie geoid undulation using (S) 
requires that CI.M l)c that value including the mass of the atmosphere. 

In the othi r (>hases i>f tlie geoid computations we must carry out the compu- 
tations aliei' the mass of the atmosphere h;is l)cen removetl from c»)nsiilerat Ion, 
nnd then add am net effects back in. To solve this pr«>blcm Moritz has shown 
that the mass <»t tlie atmosphere can be condensed «>nto the reference ellipsoid, so 
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that tho mass of the laUer will Ik- t-(|ual t»» the mass of tho solid i-arth and oceans 
plus the atinusplicre. In this ease equation (2), (12), and (17) are valid as they 
stand. To sec this more clearly, we write the Rravlty |H>tentlal at a jioint In the 
atmosphere as; 


(20) 


W \V 


r Jr r 


where M, Is the mass of the atmosphere ami M(r') Is the mass outside a sphere 
of i-adius r surrounding the earth. Kquatlon (20) is written In sliuhtly an approx- 
imated form. W Is tlic Kravity potential due only to the solid earth plus oceans. 
The normal potential Is now delined (appn>xlmately) as; 


( 21 ) 


U = -t 


r 


where \f is the norn .< jA.ientlal due to a reforenee ellipsoid that has a mass 
c((ual to that of the solid earth jjIus oceans. The disturbing potential, T, Is: 


( 22 ) 


= \v - u = \v°- ir’- G J^- 


M(r') 

r'‘ 


dr' 


(23) 


T 


= T-’ - G f” -^^dr' 
J, r' 


where. In terms of iwtential coefficients: 


(24) 


T® = 


GM 






(C, cosmX 


£ a 


:0 



Sin mX) P. (sinij). 


The scc(nid term on the right hand side of (23) can be evaluated numerically. For 




the ease of r referrinj; to the surface of the ellipsoid, the sec«>nd term has a value 
of ui)ouf 0.006 ktjal meter which wouM be alj<»ut O.G cm in ijcoid heit^ht. Thus this 
term is ncKliRlbla (’onse<iuen»ly, we can use (2H dividing by y to ol)tain N (from 
|K)(en{ial coefficients). The resultant e-Kiiression is (2). 

The compuUitior. of the Rravi’y anomaly in the case that the atmosphere is con- 
densed onto the e!li;>soid can lx? exjiressed as: 


(25) Ag = + 6g. 

where: 

. .would be the gravity anomaly in a system where all atmosphere nms.s 
is condensed info 'he ellipsoid; 

Ag would be the usual anomaly where the atmosphei’C is not condensed 

into the ellipsoid. '■* 

fig* is a correction that can be computed knowing the elevation of the 

])oint and a model for the atmosphere. 


In terms of potential coefficients, Ag^ is given by (12), or by degree in terms 
of (17). Consefjuently, the computation of in (10) or Ag, in (15) is not af- 
fected by the condensation of the atmosphere mass to the ellipsoid. 

The anomaly to be used in the Stokes' integration based on terrestial data, must 
be one to which the atmosphere correction term has been applied. Specifically, we 
must use in (10) and (. ;: 

(26) Ap ^ Ag - fig* 

where: 

will be the usually given free-air anomaly and; 

will be given based on the mean height of the (in our case) 

1®.\ 1^ l)lwk. 

Values of fig* car. be taken as -fiig where fig is tabulated for the Cicodetlc 
Heference System 1067 (l.'G, I'JTl). If we assume that fig* is constant within the 
integration cap used in the .‘blokes' integration, the effect ol the atmosphere «»n the 
gcoid undulation is: 




(27) 6N* « ^ S(U)cKt 

4ri G 

Oc 


Since (27) Is evaluated in a circular cap whose radius is ^o* we can write (27) 
in the form: 


(28) 


6N* 


“R6ir. |V 
2G 


S(C) sini;dC) 


Equation (28) can be written in the form: 


(29) 6N* “ t(tb) 

G 

where: 

(30) ?(C) = i fs(U Sint dti 


which is evaluated and tabulated in Lambert and Darlintr (1936), Using these values 
of J(t:)and values of from the Geodetic Reference System 1967; the following 
6N* values were computed. 


Table One. Effect of Atmo jphere on Ccoid Undulation Computations (meters) 


Mean Elevation of Caj) (meters) 



0 

100 

HSEHi 

300 

400 

Hi 

0. 53 

0. 56 

0.55 

0.55 

0. 54 

10' 

1.17 

1.16 

1.14 

1. IL 

1.12 

15” 

1.75 

1.73 

1.71 

1 . 69 


20” 

H5H 

mH 

2.*>1 

2. 18 

.M6 

ilSSHi 

2.C7 

2.6-1 

2.61 

WBEm 

mmm 

11 

2.97 

2.93 

2.90 

2.87 

2. 83 i 






































The corrcftions in Tabic ()nc shouI<l be added to (he undulation component com- 
puted from (10) ;md (18) wljen .\{i is taken simply as the observed free air nno- 
mally, .Tii. 

From Table (Jnc we sec that the correction is quite sensitvc to the cap size, 
l»ut insensitive to the elevation within the ranne tested. 

In tl)c actual (;coid computations to be made for this re|K>rt, 70 *a, of the gravity 
material is in ocean areas where the elevation would l>e zero. The rcmaininj{ 
areas would be in the eastem and central I'nited States, which have a small mean 
elevation. (!on.sc(|Ucut!y, for this rept)rt we will use the 6N» value for a mean 
elevation of 0 meters. A more accurate pr«)ccdure w«>uld require knowledge of 
the 1" X 1" mean elevations. lA»r this re|Mirt this Is not necessary, IjuI for liiRh 
elevation continental artas, such accurate computation for various cap sizes 
needs to t>e done in the tutu re. 


A, Numerical Integra t ion of Stokes' K(|uation 

The evaluation of (10) or (18) is carried out by summation in which an 
averufre S(V«) value is needed. We write (18) (for exaini)lc) as? 


N. y A?S(v)Aa 


where: 

Aa is the; area of the 1® x 1° block and: 

(32) S(C) - JJ S(i')fia 

Aa 




■ ! 
I 

\ 


i 


* 

I 


t ‘ 
« 


f 



i 


In practice (32) is evaluated by computing one, or several values of S(t;) from 
the comjjutatiiin |)oint to points in the block in which is (iiven, and meaning 
the result. The number of values meaned will dc)K'nd on the size of v and the 
accuracy desii-ed. Since .S(C') chanjics rai»idly for small v* values, more points 
are required for :ui accurau mean wlu-n v ia ainall than when p is lar)ce. 
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To Investlpate a proper Infci^ratlim proccHurc, \vc equate the undulation (in 
■equation form) obtained by Method A and Metliod H. The followin); equality 
Ehoidd then hold: 

,33, I (4-5 ,l-<) ) I ,C-j.cosmA 

Ob i « 

* S, sininA)'P (siif?), 

h kM 


Usi^ the GEM 6 potential coeficlents to A=lfi, the rii;ht hand eido of (33) 
can riRorously be c(»mputcd yicldinj; a "true" value. 'I'hc left hand side of (33) 
can be e\’aluatcd by subdivision schemes until one is found tliat yields results 
consistent with the "true" value. I'our subdivision schemes were tested: 


Stokes' function is computed only for the center 
point of each anomaly block (the undulation com- 
putation |)oint is at the corner of a l°x 1® block); 

Inside a s|)hcrical cap of tI)-2® around the compu- 
tation iKiint every l®x 1® block is divided into four 
equal (in terms of latitudc/longitude increments) 
blocks for which the Stokes' function is evaluated 
at tl:e center point and meaned. CXitside 
the Stokes' function is evaluated based only on the 
center jx>int of the l®x 1° block. 

Inside a spherical cap of 1,1)- 2° a 16 sub-block 
system is used, while beyond 2° a -1 sub-block 
system is implemented. 

I Sulxlivision Four (G-1-16-4-1): The followin;; sub-block system is used here: 

I 

0^2®; C-1 sub-blocks 
2®< ^l) < 5®; 16 sub-blocks 
5®< ;!) s 10®; •! sub-blocks 
10°< 0 20° ; 1 sub-block 

Usinj)^ these various sub-division scliemes, (lie loft hand side of (33) was evaluated 

i 

I • _ii- 

I 


i Subdivision One (1): 

i 

i 

I Subdivision Two (4-1) : 

% 

i 

I 

I 

t 

I 

t 

I 

t Sulxlivision Three (16-4): 

i 

I 

i 



T 


1 


1 

- T I 

i 



r'" 

"1 — 



1 

1 i 




1 


Willi the i-csulls shown in Table Two, aloni' with the "true" value coininited from 
the riiiht hand side of (3:)), for four imints in the calibration test area. 

Ta b le Two . Influence on rndulation (*om|»utat ion of Various Subdivision Schemes 
Used in the Kvaluation of the Stokes' ^'ornulla. (meters) 


Test Point (o, X) 


Sulxliv.. Model 

1 in', 277'' ) 

2(40®, 207*') 

:t(20°, 277') 

4,20' 297®) 

1 

-IS. 1 

-20. 8 

4.0 

-52. 1 

2 

-IH. 1 

-20. 1) 

4. 1 

-52. 5 

• a 

-18.4 

-27.0 

t. 1 

-52.8 

4 

-18.5 

-27.0 

4.2 

-52.9 

"true" 

-IS. 5 

-27. 1 

4.2 

-53.0 


The best agreement with the "true" value exists for subdivision model four, 
which is thus chosen for use in the actual geold (•omputalions in the calibration 
area. Since the actual integration error will depend on the magnitude of the ano- 
malies in the Stokes' kcrnal, \vc would exjiccl Method A would be somewhat less 
sensitive to integration errors tbanMethiKl H, since in A an anomaly differem-e 
is used while in Method B, the actual anomaly is used. The exact sensitivity 
will depend on the magnitude and smoothness of Ag, used in A (or in other 
words on the value of f-,*). Conscxpiently, w’e recommenit the above subdivision 
model for use in either Methotl A or Method B of undulation computation. 


5. .\ccuracv Analvsis 


5.1 Introduction 

In this section we will try to estimate the accuracy of the geoids computed by 
the methiKls described in section two. In such an analysis we will try to estimate 
the optimum truncation angle m. Several of the error terms have been discussed 
previously (Haj>p, I97:t) but they w ill lx* rediscussed here for the specific area of 
com))Utation. 

There are basically two data sixirce errors. These arise from the errors 
in the potential coefficients, and the errors in the gravity anomalies. The varl- 
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ance-co;-ariance matrix of the potential coefficients Is desljpiated Ea,, and that 
of the anomalies as Zijj. We will assume that these matrices are diaRonal. 
This assumption is true for the gravity anomalies since they are estimated in- 
dejM/ndently, but not exactly true for the potential coefficients. 


5.2 Analysis for Method A 


Method A Is rcpix*sented by (1). The first two terms on the right hand side 
of (1) represent errors of commission, while Nj represents an error of ommis- 
sion. For a single undulation error varltmce we write; 


(34) 


of + 




where: 


a, is associated with the commission errors. 

e is associated w’ith the ommission errors. 


There should be no confusion with this o. and that used previously to dcslfpiato 
the integration cap. To compute O; we write: 

(35) N# = Ni + Ns 


Considering (2) and (10), (35) can be written in matrix form as; 


(36) 


Ne = B' C + - ^.) 


where the elements of H are the coefficients of the potential coefficients 
(designated C in (36) ). The elements of A are taken from (10). We can repro 
sent the computation of agj as: 


(37) 


Ajg. ^ lie 
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T 


T 


where the elements of F}. arc found as the coefficients of the potential 
coefficients in tl2(. Insertiru; (37) into (3G) we have; 


K, » DC ^ 


Applying the propagation of error formulas we have: 


Os *1^^. n * A * 


- 23" T^, D: A. 


The evaluation of this eejuation is complex and will not Ijc attempted in this report, 
as a simpler formula will tie found when using Method H for the undulation cnmjni- 
tation. 

The specific evaluation of N, can not be done as we do not know the high 
degre<‘ poienMal coefficien's. At licst we can find a global average effect by 
writing (Rapp, TJ73): 


o| = 


■f' 


where Cffig) are anomaly degiee variances and s is (Tschernlng and Rapp, 

1971): 0.999017. Although the summation is taken to *■, in practice we would 
take the summation to about 200 as wc are using l' data as the smallest terres- 
trial data block. For e-10°, we get for about 1.1 meters and for t 20^’, 

.about 0. 7 meters. 

Cither error sources in the computation of point undulation include (for the com- 
putation in this paper) the fact that smaller anomalies than l®x 1"' anomalies may 
be needed, and that the reference set of constants may yield an equatorial gravity 
different from that implied by the best set of constants. T! ;ac effects have lieen dis- 
cussed in Rapp (1973i. For our jjurjvises wc will assume that I'^x data is (he only 
data used; (there appears to be about tO. 1 m of more detailed information in ano- 
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maly data in blocks smaller than l®x 1®); and we will assume our equatorial 
gravity’ is tnie; if it is not, errors up to 2. Cm(ti =» 20*^ ) for a 1 mpal error in 
equatorial p*avity can be expected, exclusive of No considerations (see sec, 9). 


f). 3 Analysis for Method B 


As in Method A, we have for Method B errors of commission and errors 
of ommission. The commission errors arise in (15) and (17) due to potential 
coefficient errors, and in (18) due to anomaly errors. The ommission error 
is caused by the fact that tlic summation in (15) is ttUicn to instead of*. 

Th 8 latter error is the same as ^Iveu in ( ID) and the discussion given with rc- 
sp'Xit to (‘10) is equally valid here. 

To carry out this aiialysis we express tlx; computational procedure in the 
form: 


(41) 


Nc = ifcjC + A'^ 


where: 

......... A, B, and C are as before and Q is a diagonal matrix whose 

elements are composed of values of Cj;(t!>o)> Error propagation thru* (41) 
yields for Method B: 


(42) 


a, = I/g r^, -t A^ZAgA 


Comparing (42) with the corresponding result (39), for method A indicates 
that the error analysis for method B is considerably simpler than for method 
A. The actual implementation of (42) can be done through minor modifica- 
tions oI programs used in the evaluation of (15) and (18) since and Z<lg 
are to be regarded as diagonal matrices. The accuracy of the geoid we will 
compute will be evaluated using (42) and (40) in (34). 
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G. Data Accuracy 


f 

C.l Gravity Anomaly Accuracy 




The gravity data to be used in these computations will Ik* Px 1® mean ano- 
malies. Each of these anomalies tliat has l)cen estimated has an assigned stan- 
dard deviation. In some cases an anomaly will be needed, within a cap, for 
which no estimate of that anomaly is available. In this case we will assume that 
the anomaly is zero and that its standard deviation is i 30 mgals, which is the 
root moan square variation of l®x 1® mean free-air gravity anomalies. 


The 1*^ l®data distribution witnin and around the calibration area is shown 
in Figure One. In this Figure One, zero (0) indicates those blocks for which an 
anomaly estimate is available and an * indicates tluit no anomaly exists. The 
outer borders of a 10® cap and a 20® cap are also shown. 


C. 2 Potential Ct>cfficicnt Accuracy 

The formal sLindard deviations for the GEM G potential coefficients w'cro 
provided to us l)y Frank Lcrch. These formal statistics are considered opti- 
mistic and we therefore used two approaches in estimating realistic standard 
deviations for the potential coefficients. 

In the first approach the formal standard deviations were multiplied by a 
scale factor of 3.4 as suggested by Lcrch, ct. als. (1974). An average per- 
centage error, by degree, can be uefined as: 


(43) 


" “uTiJ(^,s)r 


where: 

RMS (nie,,) ... is the root mean square value of tlie standard deviation of 
the coefficients of degree t’, and, 

RMS(C,S) . . .is the root mean square value of the cofficlents of the GEMG 
solution of degree f. 


Values of (43) are plotted in Figure Two where it is seen that (ho percentage 
error increases so tint it is about M0% at £ IG, and even higher for the c(K*f- 
ficlents above 1(J to 22. 
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Figure 1>'o, 






Wo next conskter the eominisskm error?? clue to the potential coefftcion'B. 
Again values have heen eomputecl at nine (xiints and averaged to obtain a re- 
presentative nuniljor. .Speeifk-olly, we have made coin|>utuliuns for different 
£,,, values and different values for the two error models of the (lEM C so- 
lution previously diseussed in section C.2. The resulting standard deviations 
are shown in Figure I'our for error ir.cxiel one with i,,, • tf, 12, U» and 22, 
and for 1C with error mcKlel two. Values are also given in Table Three. 

In eomjiuting the values for the 22 solution it is necessary to adopt 

standard deviations for the (uitential coefficients not present in the GK.M C set 
as we have in essence assunu*d these coefficients to be zero. Specifically, we 
have taken for their standard deviation, the i*oot mean sc]uare coefficient ^'arl- 
ation Imjdled b>' the following anomaly degree variance model given in 
Tnchcrning and Itapp (1U74, p. 20): 


(41) 


_ A^_/- IJ 

• U) 


with A - 425.28 mgal*, and B -24. 

Table 'I'hrc e. Undulation Standard Deviations IXie to Potential Coefficient 
Errors. (meters) 
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8 

* 

12 


IG 

e 00 

Error 

1 

— 

2 

Error 

1 

MocTeT 

2 

Error 

1 

Model 

2 

Error Model 
1 

0 

1.3 

3.7 

2.0 

5.0 

2.8 

5.7 

3.7 

2. 5 

1.0 

2.8 

1.3 

3.4 

1.5 

3.C 

1.0 

5.0 

0.7 

l.U 

0.7 

2.0 

0.8 

2.0 

1,3 

7.5 

0.5 

1.2 

O.G 

1.4 

1.1 

1.8 

1.7 

10.0 

0.4 

0.9 

0.7 

1.0 

1.2 

2.1 

1.4 

12.5 

0.4 

1.0 

0.8 

1.0 

1.1 

2.1 

1.2 

15.0 

0.4 

1.2 

0.8 

1.9 
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2.0 

1.1 

17. 5 

0. 1 

1.3 

O.G 

1.7 

0.7 

1.7 

1.0 

20.0 

0.4 

1.2 

0.5 

1.4 

0.7 

1.5 

O.H 

22. 5 

0.4 

1.1 

0.4 

1.2 

0.0 

1.3 

0.7 

25.0 

0.3 

0.1) 

0.1 

1.0 

0.5 

1.1 

0.0 

27.5 

0.3 

O.H 

0.3 

0.9 

0.4 

0.9 

0.5 

ao., 0 . 



.0,0 

. _!kl 

_Q,.7 

0.3 


J 





The (liffcrcnoes in the results from the two error motlcls lave a maximum 
at 0 = 0® and tleereases as 0 inereases. 

Th(? last error sourc’e considered here is that due to the ne^jlcct of the hijjh- 
er dcy;ree potential coefficients. This error is ^iven by wlu’re the summa- 
tion is taken to 200 and (-11) was used for the anomaly de^^ree variance model. 
These results are tjjven in Table Four and Fit^ure Five for various and 
values. 


Table F?^ ir. Undulation Kr*'or Uaused by Potential Uoeffieient Tiiincatlon at 
Specified (meters) 



• 




8 

12 

IG 

22 

6 

7.C 

5. 5 

4.4 

3.3 

2.5 

3.5 

1.8 

1.2 

1.0 

5.0 

2.1 

1.9 

1.8 

1.4 

7.5 

2.G 

2.3 

1.7 

0.9 


2.9 

2. 1 

1.2 

0.8 

12.5 

2.9 

1.6 

1.0 

0.8 

15.0 

2.5 

1. 1 

1.0 

O.G 

17.5 

2.0 

1.1 

0.9 

0.5 

20.0 

1.5 

1.0 

0.7 

0.5 

22.5 

1.2 

1.0 

0.5 

0.4 

2;i.O 

1.0 

0.8 

0.5 

0.3 

27.5 

0.8 

0.5 

0.4 

0.3 

30.0 

0.7 

0.1 

0.3 

0.2 


We can now quadratienlly add the three error sourees topethcr to obtain the 
final ix'prt\scnta( ivc undulation standard deviation in the calibration area. The 
resultant st:uuLird deviation'; arc j)lotted in Fij^ure Six. 

In considerinj; Fii^jure Six tor the S, 12, 10, and 22 solutions with |X)tential 
coefficient crr«)r model one we see that the error at first decreases as v in- 
creases, then it increases somewhat and finally decreases. At certain R values 
such no 7. 5® for -f.,, 12 we would e>qx*et a Iart;er error than at smaller t>r 

larger C' vahies. This would indicate that one should choose an ojdimum c va- 
lue based on the of the potential eoeffieien’ field being used. lUised on 
Figure Six we feel a i., -O' cap is reasonable. 


Pace 
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Figure .Six . nepresentaHvf Tola’ I’ndulation .Standard Pcviatlon in the 
Calilji-ation Area. 

with PotonMal (’<K*fficic‘nt M»>dcl one for £t»» 

12, IG, and 22 

wi*h PofenHal Coefficient Mixlel Two for L»tt ■ IG 


The st.andard deviation of the undulation when the pc.ssiinistic error model 
for the potential coeffieients is used is shown in Figure Six for the case of f.,, 
= IG, Although the errors are larger, the difference decrea.ses as c increases 
so that at v 20” and 16, the total standard deviation of the undulation us 

ing the oj)timistic potential ct»efficient error model is il. 5m while it is s2. 1 m 
when using tlie pessimistic error model. 





For the final undulation computations the GEM f> solution was used with 
= 10. This was done so tliat the error analysis could be done without 
making assumptions about the coefficients between degree 17 and 22 not es- 
timated in the GEM 0 solution. It will be shown that no significant undula- 
tion difference occurs when using a ^16 or 22 in the case of the GEM C 
set. 


8. The Geoid in the Calibration Area 

The previous discussions have dealt with the method for detailed undu- 
lation computation and the method of accuracy analysis. \Vc now turn to the 
computation of the detailed geoid and its accuracy in the calibration area us- 
ing the GEM C poten'ial coefficients and 1° gravity anomalies given to a 
mgtil on a tape dated July li)75. 

W’c first define a set of constants identical to that used by Marsh and 
Vincent (1973) so that our undulations may be compared. These constants are: 


o;= 7.2921151467 x lO'^^rad/s 
a = 6376142m 
f=- 1/296.255 

GM= 3.966009 X 10''m=*/s® 

Wo 6263667.52 krai m 
y, = 97S032. 14 mgal 


The GM includes the mass of the atmosnhere*nnd thus the effect of the at- 
mosphere is included in the geoid potential and equatorial gravity, y, . The 
normal gravity formula corresponding to these constants is then: 


(•15) = y.(l ^ 0.0053024269 sin^O - 0.0000059 8in‘"2o) 


Now the anomalies on the July 1975, 1® tape were given with respect to the 
gravity formula of the Geodetic Hefercnce System 1967 and thus need to be con- 
verted to be given with rcs|)ect to the new constants. We have: 


♦ If it did not, then the Nn term from equation (51) would have to be evaluaied. 


(IG) 


^BiS7 * >(5 7 “ Vr** 




or numerically: 


(17) ig„, = - 0.29 - O.OGsin^o nii;j(l 


Although this correction is small, it is systematic and its neglect could cause 
errors on the order of 0.8 meters in the computed undulation. 

For the final results several different computations were performed us- 
ing the GEM C potential coefficients and the 1®.\ 1® mean anomalies converted 
to a gravity formula consistent with the adopted constants. We first-computed 
the geoid undulations using Method B with the GEM 6 coefficients truncated 
at degree IG and with a cap size of 20®. These undulations arc given in Table 
Five and in contour form in Figure Seven. These undulations include an at- 
mospheric correction of 2.2C. m obtained from Table One. The undulation stan- 
dard deviations for the calibration area when using the potential coefficient 
error model one is shown in Table .Six and when using the potential coefficient 
error model two in Table Seven. In the first case the standard deviations have 
a maximum of r2.0 m with a minimum standard deviation of tl. 1 m. The cor- 
responding values when usmg the second error model are r2.-l m and i-l.Sm. 
These error estimates are all given with respect to the adopted set of constants, 
and exclude the error contribution due to the neglect i>f detailed vcravity informa- 
tion in blocks ix'low l®x 1® in size. The accuracy of these standard deviations 
depends on hou well the accuracy models for the anomalies, potential coefficients 
and truncation effects liavc been handled. We believe that the latter two effects 
have Ix^en handled reasonably. However, the anomaly error contribution has 
been based on the anomaly standard deviations that may be optimistic. 

We next computed the geoid undulations using .Method A described in sec- 
tion 2. 1 with the identical data as used in .Method B. A comparison of the n sul- 
tant geoids showed a maximum descrc|<*ncy of 0.2 m with other statistics on the 
eomjiarison given in Table Eight. Wc ciinclude that either method consistently 
apjjlied will give the same undulation at the :U. 1 meter level. 

\ comjHitation was m.atle using method B when a truncation angle of lO® was 
used in.stcad of 20®. The largest difli-i enci' found was .9. f* meters with the 
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square rwt of the variance of (he difference beinij <-1. 55 meters. Hased on 
the stantlard doviati«)ns shonn in l-’ii,r>iro Six, we would ex|)ec( the above dif- 
ference to l>c when usinj; the optimistic potential coefficient error 

mo<lcl and tl. Sin when ui^inti the pessimistic model. (The values are arrived 
at by computin;; - a-V - o^r'i. These exi>ected values are very com- 

patible with the result actually found. 

In the computation described in the above pai*actraph, the atmospheric cor- 
rection was riRorousIy applied with mean difference between the two results l>e- 
inj{ -0.24 meters. A test comj'Utation was made when the atmosphere was not 
considered in the anomaly data within the 10’ and 20® caps. In this case the 
rcsultinit mean undulation difference was found to be -1..14 meters which is a 
considerable increase from the -0.21 meters when tlie atmosphere was proper- 
ly treated. These results indicate the practical value and need of the atmo- 
spheric correction to {jravity anomalies when the undulations of the ijeold are 
being con.puted. 

Another conqxitation that was made was with the GEM fi p<)(ential coeffi- 
cients tal'.cn to I = 22 with a u 20®. The resulting undulations were conq»arod 
to the values computed using the coefficients to I IG. The differences l)c- 
tween the undulations was a maximum of 0.3 m with a negligible mean differ- 
ence imd variance as is shown in Table Eight. 


Table Eight . Comparison of Various Undulation Computations (meters) 


1 

A 

B 

C' 

D 

E 

Mean Diff 

.03 

-0.21 

.03 

3.87 

•1. 15 

a 

i.09 

1 1 . 55 

±. IH 

t2.59 

±3. 52 

Max (A Diff 

0.2 

2.9 

0.3 

12.0 

10. 1 

Max W Diff 

-0.2 

-3.6 

-0.3 

-0.6 

-5. 3 


A: Difference between Method A minus Method B (*,,, 16, 0 20®); 

B: Difference bemeen Method B(C 20®) minus Method B (W=- 10®); 

C: Difference between Method Bf*,,, 22) minus Method B(i'„, 16); 

D: Diffci*encc behveen Method B (t. 20®)minus solution of Vincent-Marsh; 

E: Difference between Mcthi'd B (t 10®)minus solution of Vincent-Marsh, 

The final comparison that was carried out was the comparison of the cali- 
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bratlon of this rO|K(rl to that ;jpoI( 1 of Vincent and Marsh (1973), which 
was provided tt> us by Jim Marsh at Fx P corner points. The results of the 
comparison are shown in Table Eij;ht. Here we see (hat there is a strong 
systematic difference in (he undulation of 3. S7 m, about 1. ? meters of which 
is caused by the nc' 4 lcct liy Vincent and Marsh of the influence of the atmos- 
phere. The square root of the variance of the undulation differences is *2. 
meters with the largc-st discrcijancy being 12.0 m. (Again |>art of this latter 
diserepancy is due to their neglect ol the atmosphere). The main cause »)f 
the non-systematic differences between the undulations is probably the Px 1® 
gravity data used in the computations. 


0. The Zero-t'irder fnchilaijo n 

The zenj-ordcr urdulation of the geoid has been discussed by Ilciskanen 
and Moritz (1907, p. 102) for the case of an undulation computation using Stokes 
cc{uation in a global integration process. .Specifically we have: 


(-18) 




k 6 m 

2 G H 


2G 


where k 6 m is the differeive berween a (laie value of the geocentric gravitational 
constant of the earth plus the atmosphere, and that adopted for the reference field. 
In addition jgj is the mean gravity anomah, of those anomalies (after the atmos- 
pheric cerroction has been applied) referred to the adopted eoustnuis. Thus: 

(•19) r <.ig-6g.)da 


This Jigj must also be subtracted from the given in (2G) to assure that 
the anomalies used in Stokes' equation liave a global average equal to zero. In 
(his case an additional term will appear from (10) and (lb) equal to: 


( 50 ) 


•ItG 






Qt-4 




r 

^ S;w) da 
(Tc 


G 
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Equation (50) is now added »o (-IS) to obtain the imxlified N, term that applies 
for both methods of computation of ^^eoid undulations described in this report. 
We have: 


(51, 


2 G It G 



If is IbO® we have » N;. If kfiM and arc zero Nc is zero. If 
R6 m b 0, and -g. * Imgal, N, is -4. Cm for Cj * 10® and -5.8m for C-, -20®. 
The value of <lg; can be determined from (49) if we have a global gravity field 
or from: 


(52) Ago y - y 

where y’ is the true equatorial gravity based on the mass of the earth plus the 
atmosph and y is the corresponding value adopted for use in the computa- 
tions. 

Thus, we have derived a new N- term that is cap size dependent. This 
No. If it can be determined, must be added to the undulations obtained from 
equation il) or equation (Mi to obtain the true undulation with respect »o the 
adopted constants, if N, is set to zero, the coinixitcd undulations will refer to 
a set of unknov.-n constants such that k 6 M and Agj are zero. 


10. Summary 

This paper documents two procedures that can be used for the computa- 
tion of geoid undulations combining po*cnti:il coefficient data and terrestrial 
gravity data. \Vc found that the two iucIkhIs yteld essentially the same rcsul's 
(at the r0.09ni levcli. However, the error analysis for Method H is simpler . 
than that for Method 

In developing the pi\)eedures for each method two important techniques 
must bo used. First, it was found that liie atiuo.spheric correction is slgnifi- 
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canf and cannot ho ncKloctid as Its neglect can cause errors in the computed 
undulation on the order of 2 meters with a truncation cap of 20''. Second, 
the numerical uv curat ion of S'okes' equation must be done in a precise man- 
ner or integration errors of about t meter will result. 

A preliminary error analysis was done considering he error sources 
due to potential co<?fficicnt errors, gravity anomaly error.* , and the truncation 
error caused by lal-.ing the potcnMal cocfficicn’.s to a certain maximum degree 
only. This annlvsis indicated that ccr'ain cap sizes larouiid 10"' i would give 
poorer results than more optimum sizes, such as 20®, which was selected 
for use here. 

The actual geoid undulations were cominitcd in the Geos - 3 calibration 
area using the GKM6 potential coefficients to degree 1C and the f x 1® mean 
anomalies available in .July 11*75. The results obtained showed a mean differ- 
ence of 3.87 m from the ineent -M.tr. '^h gcoicl with a difference variance of 
(2. 58 m)‘. The estimated standard deviation of the undulations computed here 
were on the order of 1 to 2 meters with respect *o the defined constants. This 
error analysis neglects the effect of using anomaly blocks of a size smaller 
than 1® X 1®, and the ellcct of ellipsoidal correction terms to Stokes' equa- 
tion. 


Although these computations have been done for the GEM 6 coeffieients 
they can easily be repeated for other eocfficient sets with a corresponding 
error analysis provided a realistic variance-covariance matrix for the poten- 
tial eoeffieic-nts is arailablc. 
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